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Abstract In the context of Markov processes, we show a new scheme to 
derive dual processes and a duality function based on a boson representation. 
This scheme is applicable to a case in which a generator is expressed by boson 
creation and annihilation operators. For some stochastic processes, duality 
relations have been known, which connect continuous time Markov processes 
with discrete state space and those with continuous state space. We clarify 
that using a generating function approach and the Doi-Peliti method, a birth- 
death process (or discrete random walk model) is naturally connected to a 
differential equation with continuous variables, which would be interpreted 
as a dual Markov process. The key point in the derivation is to use bosonic 
coherent states as a bra state, instead of a conventional projection state. As 
examples, we apply the scheme to a simple birth-coagulation process and a 
Brownian momentum process. The generator of the Brownian momentum 
process is written by elements of the SU (1,1) algebra, and using a boson 
realization of SU(1, 1) we show that the same scheme is available. 



1 Introduction 

In the context of nonequilibrium physics, simple systems of interacting par- 
ticles have been received considerable attention recently, and it has been 
known that the concept of duality is useful in studying stochastic processes 
of interacting particle systems [1]. The duality would give deep insights and 
analytical results for stochastic models in nonequilibrium, and actually, there 
are many works using duality properties, ranging from calculations of corre- 
lation functions in interacting particle systems [2-8] to studies of biological 
population models [9, 10]. 

While the duality is a useful concept, there are some problems to use it. 
For example, it has been sometimes necessary to construct a dual process with 
ad-hoc procedures. In addition, while the original process and the dual one 
are connected via a duality function, the duality function should be selected 
properly. To construct the duality function, usually one sets an ansatz for 
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the duality function, and checks whether the ansatz satisfies the dual relation 
or not. Recently, a general procedure to derive a duality function has been 
proposed [8] ; dual relations for various stochastic models have been recovered 
using the symmetries of the original process. While this general procedure 
has wide applications, it may be needed to use a heuristic way for some 
specific cases. For example, the general procedure in [8] may not be available 
for Brownian momentum (or energy) processes in boundary driven cases, and 
a boundary part in a duality function was heuristically found in [8]. 

In the present paper, we show a new scheme to obtain a dual process and 
a duality function. The scheme is based on a boson representation, and it is 
applicable to a case in which a generator of the stochastic process is expressed 
by boson creation and annihilation operators. We will see that a continuous 
time Markov process with discrete state space (e.g., simple birth-death pro- 
cesses and discrete random walk models) is a dual process of a continuous 
time Markov process with continuous state space (e.g., stochastic differential 
equations and Brownian momentum models). In the derivation of the du- 
ality function, the generating function approach and the Doi-Peliti method 
are used. It will be clarified that the dual process and the duality function 
are naturally derived by using bosonic coherent states as a bra state in the 
Doi-Peliti method, instead of a usual projection state. In addition, using a 
boson realization of SU (1, 1), it is also possible to study a duality relation for 
a stochastic model with elements of the SU(1, 1) algebra. Especially, we can 
derive a duality function for a Brownian momentum process with boundaries 
not heuristically, but deductively. 

The outline of the present paper is as follows. In section 2, we give a 
definition of duality. The new scheme to derive a dual process and a dual- 
ity function is shown in section 3. Sections 4 and 5 are applications of the 
new scheme to two examples; i.e., a simple birth-coagulation process and a 
Brownian momentum process. Section 6 gives concluding remarks. 

2 Duality 

General discussions for duality are given in [1]. In the present paper, we only 
treat the duality between birth-death processes (or discrete random walk 
models) and diffusion processes. 

Suppose that (£t)t>o and {zt)t>o are continuous time Markov processes 
on state spaces J? and J?duai> respectively. Let denotes the expectation 
given that the process (£t)t>o starts from £. The process (£t)t>o is said to be 
dual to (z t )t>o with respect to a duality function D : Q x J?duai — > M if for 
all £ G Q, z G J?duai and t > we have 

E s D(z,£ t )=E< i ™ i D(z u Z), (1) 

where Ef ual is expectation in the process (z t )t>o starting from z. 

In the following discussions and examples, the process (£t)t>o is a con- 
tinuous time Markov process denoting a birth-death process (or a discrete 
random walk model), so that £ t G N. On the other hand, the dual process 
(zt)t>o is a continuous time Markov process with continuous variables, and 
then z t el. 
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3 Derivation of duality function using boson representation 

In this section, we derive a duality function. In order to obtain it, we firstly 
explain a correspondence between a generating function approach and the 
Doi-Pcliti method (second quantization method). After that, it will be shown 
that a duality function is naturally obtained from a state vector for the 
continuous time Markov process with discrete state space. 



3.1 Generating function approach 

Some stochastic models with discrete variables are described as birth-death 
processes. A time evolution of a birth-death system obeys a master equation, 
and it is sometimes useful to treat a generating function instead of the original 
master equation [11]. For simplicity, we here treat a birth-death process with 
only one variable. The generating function G(x, t) is defined as 

oo 

G(x,t) = Y t P(n,t)x n , (2) 

n=0 

where n £ N, x € K, and P(n, t) is the probability with n particles at time 
t. The time evolution equation for G(x, t) is written as 

j t G(x,t)=L^x,^jG(x,t), (3) 

where L(x, ^) is a linear operator, which is constructed from the original 
master equation. 



3.2 Doi-Peliti method: boson representation 

The Doi-Pcliti method is a well-known approach to investigate birth-death 
systems [12-14]. In the Doi-Peliti method, bosonic creation and annihilation 
operators are used: the creation operator at and annihilation operator a 
satisfy the commutation relations 

[a,a f ] = l, [a, a] = [a 1 , a 1 ] = 0, (4) 

and each operator works on a vector in Fock space \n) as follows: 

a^\n) = \n + 1), a\n) = n\n — 1). (5) 

The vacuum state |0) is characterized by a\0) — 0. The inner product of bra 
state (ra| and ket state \n) is defined as 

(m\n) = 5 m , n n!, (6) 

where 8 m ^ n is the Kronecker delta. 
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When we define a time-dependent state \ip(t)) as 

oo 

h/>(t)) = XXM)|n>, (7) 

n=Q 

the time evolution of the state \4>(t)) is given by 

±m))=Ltf,°)W)h (8) 

which recovers the original master equation. The linear operator L(a\a) is 
obtained from the original master equation, and it is known that L(aJ ' , a) in 
© has the same form as L(x, ■4-) in 

While the Doi-Peliti method is similar with usual quantum mechanics, 
there are some differences. One of the big differences is the usage of a pro- 
jection state. In the Doi-Peliti method, the projection state 

<^E^< n i = <oie a ( 9 ) 

n=0 

is used to obtain physical quantities. For example, the average of n is given 



3.3 Connection between generating function approach and Doi-Peliti 
method 

There is a one-to-one correspondence between the generating function ap- 
proach and the Doi-Peliti method. Consider the following construction for 
ket and bra states in the Doi-Peliti method: 

\n) = x n , (m\ = J dx8(x)(J^ (•), (10) 

and interpret the creation and annihilation operators as follows: 

a) = x, a = — . (11) 
da; 

Hence, we immediately see that all properties in the Doi-Peliti method are 
recovered using x and -p. In addition, the linear operator L{a\a) in (JSj) is 
obtained by replacing x and -nr of L(x, -jM in |3]) with and a, respectively. 

Because of the correspondence between the generating function approach 
and the Doi-Peliti method, we will switch between these two notations freely 
in the following discussions. 
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3.4 Derivation of duality function 

For simplicity, a case with only one variable is discussed at first. After that, 
a result for general cases will be given. 

We consider the following time evolution equation 

±\<f>(t)) = L\m, (12) 

where \4>(t)} is given by 

oo 

|#t)> = X;P&t)|0, |O=d(at,O|0>. (13) 

Here, £ is a variable used in a continuous time Markov process with discrete 
state space (£ e N), whose probability distribution is denoted by P(£,t). 
Note that the state |£) is not restricted to the form discussed in section 3.2, 
and |£) is generated using the creation operator via a function d(a',£). 
For example, we will see the following functions d(a^,£) in sections 4 and 5: 

In section 4: d(aJ,t;) = (a*)^, 

(at) 2 « 

In section 5: d(a',£) = 



where (2n — 1)!! = (2n — l)(2n — 3) • • • 3 • 1. In both cases, the state |£) is 
expressed in terms of the creation operator a' , but it is not necessary to use 
the simple construction \n) = (a t )"|0) in section 3.2. 

As explained in section 3.2, the projection state is usually used as an 
adequate 'bra' state in the Doi-Peliti formalism. The key point to obtain a 
duality function here is the following one; instead of the projection state, we 
define a bra state {(j)(t)\ as 

/oo 
dz4>(z,t)(z\, (14) 
-OO 

where (z\ is a coherent state of o)\ 

(z\ = {0\e za , (15) 

which satisfies 



z\a) = z(z\, 



(16) 



and z is assumed to be a real variable. From the correspondence between the 
generating function approach and Doi-Peliti method, the following identities 
are easily checked: 

(,t|») = z", (17) 

eN*(sY !»>-.* (£)'{.!«>. (18) 
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The linear operator L{a),a) (i.e, L(x, 4-)) is generally written in normal 
order, i.e., all creation operators are to the left of all annihilation operators 
in products. Hence, 



(4>(t)\L[x, 



dx 



dzcf>(z, t)(z\L I x, 



dx 



dz<f>(z,t)L ( z, — ) (z\ 



dz 



L* z, 



dz 



(z,t) 



{A 



(19) 



where L (z, 4M is obtained by simply replacing x and 4= as z and 4-, respec- 
tively; L* (z, 4^) is the adjoint operator of L (z, 4-)- Therefore, we obtain 
the following identity: 



where the time development of the bra state (</>(£) | is defined as 

^<<Xz,t)| = (4>(z,t)\L. 
Combining (fP9"]l and (|2T|) . we have 

±4,(z,t)=L* (z,^) 4>(z,t). 



(20) 



(21) 



(22) 



At this stage, it is clarified that the linear operator L is a generator for the 
continuous time Markov process (zt)t>o if 4>( z it) can be considered as the 
time-dependent probability density. Hence, a continuous time Markov process 
with discrete state space, (£t)t>0j is naturally connected to the stochastic 
process with continuous variables (zt)t>o- In addition, writing ([2H)l explicitly, 
we have 

/oo 00 />oo 00 

dzJ24>(z^) p &t)d(z,0= dz^0(z,t)P(C,O)d(z,e). (23) 
-oo (=Q J-oc £=Q 

Hence, if we set the initial conditions for £ and z as a Kronecker delta function 
and a Dirac delta function respectively, the duality relation (p} is obtained. It 
is also clear that the function d gives a duality function. The above discussion 
means that if a generator L is expressed in terms of creation and annihilation 
operators, a dual process is immediately constructed and the duality function 
is given by the function d, which specifies a state £) in the Markov process 
with discrete state space (see dT^l) ). 

If a Markov process with discrete state space has many variables 
(i E {1, . . . , N}), a state |£) is defined by 



N 



N 
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U.di(aUi 



i=l 



\i=X 




(24) 



7 



where di(a], £) may be different from each other. Hence, using z = {z\, z 2 , . . . , zjy}, 
Zi G M, and £ = {£1, £2, • • • , £jv}, £,i £ N, a duality function is given as 

N 

D(z,0 = Hdi{zi,^. (25) 



4 Example 1: Simple birth-coagulation process 

As a first example, we apply the scheme in section 3 to a simple birth- 
coagulation process. The birth-coagulation process has been used widely to 
study front-propagation problems, and it has been known that a Langevin 
equation, so-called stochastic Fisher and Kolmogorov-Petrovsky-Piscounov 
(sFKPP) equation, plays an important role in the study of the front-propagation 
problems [15-18]. Recently, the sFKPP equation has been discussed even in 
a QCD context [19]. 

A duality relation for the birth-coagulation process has been used to 
study a front propagating problem in [6]. In the duality relation, the birth- 
coagulation process is connected to a Langevin equation. In [6], the dual pro- 
cess and the duality function were assumed, and explicit calculations based 
on stochastic differential equations were used to check the duality relation. 
We will show that the dual process and the duality function are recovered 
simply using our general scheme. In addition, we will derive a new duality re- 
lation for a slightly-changed stochastic process. The derivation demonstrates 
the effectiveness of the present scheme to find a new duality relation. 



4.1 Derivation of duality relation in the birth-coagulation process 
Consider the following reaction scheme for the birth-coagulation process: 

A — *■ A + A at rate 7, 

A + A^A at rate a 2 . (26) 
The master equation for the birth-coagulation process is written as 

t) =t(£ - i)p(t - 1, t) - 7 ep(£, t) 

a 2 fc^P(£, t) + a 2 ii±^P(£ + 1, t), (27) 

where £ is the number of particle A, and £ e N. The linear operator in the 
Doi-Peliti method is given by 

2 

L = 7 (a t -l)a t a+y(l-a t )a t a 2 , (28) 
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and it is easy to check that the following time evolution equation and a bra 
state \4>{t)) recover the master equation (j27jl : 

, oo 

±\m=i<w)), i^(*)>=E^,*)io (29) 

5=0 

and 

10 = (at)«|0). (30) 
Since d(aJ,£) = a duality function is 

D{z,t)=zt. (31) 
Considering the adjoint of the linear operator L in terms of z and 4-, 

L* = -± [- 7 ,(1 - Z )] + [a 2 (l - z)z] , (32) 

we see that the adjoint operator L* gives a Fokker-Planck equation. Hence, 
the dual process corresponds to the following stochastic differential equation: 

dz = -7*(1 - z)dt + ay/z{l - z)dW. (33) 

If one consider a new process via a variable transformation it(t) = 1 — z(t), 
the corresponding stochastic differential equation is 

du = 7u(l - u)dt + ay/u(l - ujdW, (34) 

and the duality function is rewritten by using the new variable u as 

D(u,0 = {l-u)£. (35) 

The above dual process and the duality function are consistent with results 
in [6]. 

4.2 New duality relation for a slightly-changed reaction scheme 

We here show that it is easy to derive a new duality relation for a slightly- 
changed reaction scheme, which has not been studied yet. 

We consider that the following reaction is added to the stochastic system 

A -> A + A + A at rate a. (36) 

Hence, a new term, a(£ — 2)P(£ — 2,t) — a£P(£,t), is added to the master 
equation (|27p . The corresponding linear operator in the Doi-Peliti method is 
a((a^) 2 — l)a^a, and then the adjoint operator L* is finally given by 

L* = ~ [- 7 z(l - z) - az(l - z 2 )] + ^ [a 2 {\ - z)z] . (37) 

Thus, we conclude that the dual process is given by the following stochastic 
differential equation: 

dz = -[72(1 - z) + az(l - z 2 )]dt + cjy/z(l - z)dW, (38) 
and the duality function is given by (I3ip . 
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5 Example 2: Brownian momentum process 

As a next example, a Brownian momentum process is studied [7,8]. In this 
case, different from example 1 in section 4, we start from a stochastic process 
with continuous state space and obtain a dual process with discrete state 
space. 

The model is defined as a stochastic process on TV-dimensional vectors 
(zi, . . . , zjy) e R N , which have to be interpreted as momenta associated with 
lattice sites {1, . . . , N}. In addition, sites 1 and N are in contact with heat 
reservoirs at temperature TL and Tr, respectively. The process is defined by 
a generator L as follows: 

JV-l 

L = L 1 +L N +Y, L i,i+i> ( 39 ) 

i=l 

with 

^f = ^f~^ (40) 

&_ f _ ZN _d_ 

" dz 2 N dz N - 

W^^-^j (/), (42) 

where / is a C°° function. Hence, the time-dependent probability density 
p(z,t) obeys the following equation: 

-p(z,t) = L*p(z,t), (43) 

N-l 

L*=L* 1+ L* N +J2Lii+i, (44) 

i=l 



LNf = T R —f-z N — f, (41) 

2 



where 



L . Nf = TR jL. f + : ±. M) , (46) 

= L w+ i. (47) 

The above equation corresponds to a stochastic process with continuous state 
space [7]. 

It has been known that the generator L is rewritten by using elements 
(operators) of the SU(1, 1) algebra [8]. The operators are defined by 

K-y. «r-i^. «?-i(^<--» + «'^)' (48) 
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and they satisfy the following commutation relations: 

[K?,K±] = ±K±, [Kr,Kf] = 2Kl (49) 
The components of the generator L is rewritten as 

= 4 (V^~ + i + KrK+ +1 2K?K? +1 + , (50) 

Li = 2T L Jff - 2K° + i (51) 

£jv = ^TrK^ — 2K% + ^. (52) 

In addition, the SU(1, 1) group admits a discrete (infinite dimensional) rep- 
resentation: 

^+16) = Q+6) I6 + 1), ^H6) = 6i6-i), x°i6) = (6 + 16). 

(53) 

Here, we reinterpret the operators L\ and Ln as follows: 

Lx = 2T L ^- - 2A? + X - = 2a\K^ - 2K* + 1 (54) 

L,v = 2TrK n -2K N + l - = 2a N+1 K N -2K% + ^, (55) 

where we interpret the constants Tl and Tr as the creation operators related 
to additional sites and N + 1; Tl = zq = aj and Tr = Zjv+i = a N+1 . 
This rcintcrprctation is justified according to the correspondence between 
the generating function approach and the Doi-Peliti method (see section 3.3). 
Using the introduction of the bosonic creation operators for sites and TV + 1 
and the SU (1,1) algebra for the other sites, it is possible to consider that the 
generator L creates a stochastic process with discrete state space. We define, 
for £ G f2, i, j G {0, . . . , N + 1}, the configuration £ lJ to be the configuration 
obtained from £ by removing one particle at i and adding one particle at j. 
Hence, the linear operator L is interpreted as 

14(0 = 

2&M£ 1,0 ) - + 26(26 + 1)^ 1 ' 2 ) - VK6)] 

JV-l 

+ £ (26(26_! + 1M6' 4 - 1 ) - V(01 + 26(26+1 + iM6' m ) - ^(6)]) 

»=2 

+ 2fcv(2&v-i + l)^"" 1 ) - + 26v[V(£ JV ' JV+1 ) - V(6)L (56) 

where ip : Q — ► R is an arbitrary function of the finite particle configurations. 
This process is considered as a discrete random walk model with absorbing 
sites and N + 1 [7]. Hence, a stochastic process with continuous variables 
(the Brownian momentum process) is naturally connected to a stochastic 
process with discrete variables (discrete random walk model). 
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Next, we obtain the duality function. Using (l53|) iteratively, it is easy to 
confirm that a state for i £ {1 , . . . , N} is given by 



l&> = 



2- 



(2& - 1)!! 



(57) 



where |0)£ is the vacuum state in the SU(1, 1) representation. Here, we note 
that there are boson representations for the operators in SU(1, 1) [20]: 



o(«!) 2 , K: 



\(^)\ K° = ^4+4^). (58) 



(See also Hence, it is possible to rewrite state by using boson 

creation operators instead of Kf . On the other hand, for site and iV+1, it is 
necessary to indicate a state |£o)o (|£iv+i)iv+i) by using the creation operator 
°l ( a jv+i) instead of K$ (K^ +1 ); i.e., |£ )o = (4) e °l°)o and |£iv+i)jv+i = 
( a ]v+i) €jv+1 |0}at+i- We therefore obtain 



N 



io = i£o)o ®ie»>» ® ieiv+i)iv+i 



((4) Co |0)o 



JV 



(26 - 1)! 



(al) 2e< |0), 



® ((aWi) C " +1 |0>JV+i) 



This means 



f(o|)«« fori = 0,JV+l, 



otherwise, 



(59) 



(60) 



and then we immediately have a duality function as (note that zq = Tl and 

zn+i = Tr) 



N 



Mi 



(61) 



which is consistent with the results of [7, 



6 Concluding remarks 

In the present paper, a new scheme to obtain a dual process and duality 
function was given; the scheme is applicable when a generator is expressed in 
terms of boson representation. Using the new scheme, it is possible to find the 
duality function not heuristically, but deductively using a state expression 
for a continuous time Markov process with discrete state space. The duality 
connects the Markov process to the dual process with continuous state space. 
We applied the scheme to two examples, and adequately recovered results in 
previous works, and derived a new duality relation for a slightly-changed 
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stochastic process. We here note that the same scheme is available to obtain 
a duality function in a Brownian energy model in [8] , and we confirmed that 
a suitable duality function is actually obtained. 

Our analysis in the present paper was limited to simple cases in which 
a representation of the x-representation is changed to the n-representation. 
In [8], the role of symmetries was considered, which then leads to self-duality 
for discrete processes. The key idea in our scheme is a usage of coherent 
states as a bra state, instead of a conventional projection state. Although 
we here limited ourselves to cases with boson representation, the similar 
idea will be also applicable to other cases. For example, it will be possible 
to apply the similar discussions to a case with a generator with the SU(2) 
algebra (or the quantum algebra U q [SU(2)]), which describes symmetric (or 
asymmetric) simple exclusion processes [5], Such extensions are out of the 
scope of the present paper; these issues are currently under investigations 
and will be published in future. 

Finally, we comment on the applicability of the scheme. The scheme would 
basically connect a discrete representation and continuous one, and give a 
duality function between them. However, there may be no guarantee that a 
corresponding dual 'process' adequately describes a stochastic process; e.g., 
it is necessary for Markov processes with discrete state space to have the 
generators with some specific forms (see (|28)l and (j56|) ). In our scheme, it is 
at least possible to derive a differential equation with continuous variables, 
which may not be interpreted as a stochastic process. The derived differen- 
tial equation and the original Markov process is adequately connected via 
a (duality) function, which is simply given from a state expression for the 
Markov process, as discussed in the present paper. Such 'duality' between 
a stochastic process and a deterministic differential equation would be valu- 
able for studies of noncquilibrium physics. In addition, the present formalism 
based on the generating function and Doi-Peliti method would be tractable 
for physicists, and then it will help to seek new duality relations for stochastic 
processes. 
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